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The Lindblad generators of the master equation define which kind of decoherence happens in an 
open quantum system. We are working with a two qubit system and choose the generators to be 
projection operators on the eigenstates of the system and unitary bilocal rotations of them. The 
resulting decoherence modes are studied in detail. Besides the general solutions we investigate the 
special case of maximally entangled states - the Bell singlet states. The results are depicted in the 
so-called spin geometry picture which allows to illustrate the evolution of the (nonlocal) correlations 
stored in a certain state. The question for which conditions the path traced out in the geometric 
picture depends only on the relative angle between the bilocal rotations is addressed. 
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I. INTRODUCTION 

The theory of open quantum systems plays a major role in many applications of quantum physics since 
perfect isolation of a quantum system is never possible. Because the environmental degrees of freedom 
are not accessible the dynamics of open quantum systems are described by effective dynamics: the 
quantum master equation [l| . The notion of decoherence is introduced which describes the loss of quantum 
coherences in a system coupled to an external environment. This concept allows to understand the 
transition from a quantum to a classical world [3] . 

Geometric pictures of the quantum mechanical state space of a system are attracting and illustrative and 
provide deeper insight into unsolved problems. Therefore it is tempting to find geometric pictures for the 
whole state space of density matrices which is nearly impossible because even for the simplest case of a 
qubit the whole state space has 4 real dimensions. 

By restricting to several properties one can find appealing pictures and visualizations of the state space 
(see for instance Q). E.g., for pure two qubit states a Hopf map can be found which provides an 
entanglement sensitive stratification of the state space Q. The so-called spin geometric picture allows 
for another visualization of the two qubit system and was introduced and discussed by the Horodeckis 
0,0], VoUbrecht and Werner Q and Bertlmann, Narnhofer and Thirring 

The paper is organized in the following way. The next subsections give a brief introduction to the theory 
of decoherence with special emphasis on the theoretical formulation used later on in this paper and an 
introduction to the spin geometry picture where also the role of the singular value decomposition is 
emphasized. In Sect. lUwe introduce the decoherence modes under consideration which arise by local 
unitary rotations of the projection operators on the eigenstates of the system. We present the general 
solution of the time evolution of the local parameters m and n and the correlation matrix c for three 
different types of decoherence, which we call decoherence modes. These general solutions are illuminated 
by considering special initial conditions: the Bell singlet state (Sect. IIIII) . The evolution of the states 
under the investigated decoherence modes is visualized graphically within the spin geometry picture. In 
the last section IIVI we show which restrictions on the initial correlation matrix c of a state with local 
parameters equal to zero have to be satisfied such that the geometric pictures of mode B and C are equal. 
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A. Decoherence 



A quantum system S is coupled to the environment E and the closed total system S + E is governed by 
unitary evolution given by the Hamilton operator Hs+sit) = H(^l + 1(^ He + Hj, where H and He are 
the free Hamiltonians of the system and the environment, respectively, and Hj denotes the interaction 
Hamiltonian. Under several assumptions (see, e.g., [l|) the reduced non-unitary dynamics of the open 
system S are given by the master equation 

^^p{t)^-^[H,pit)]-Vip{t)), (1) 

with the dissipator [1, 

^(P) = ^ E i^l^f'P + P^l^^ - '^^^pA) ■ (2) 

k 

The general structure of the dissipator can be simplified by assuming the a priori arbitrary Lindblad 
generators A], to be projection operators A}^ — ^/\kPk, with P| = Pk (see Ref. [HI), which gives a 
simplified structure of the dissipator 

V{p) = \Y.^>' (^''P + PP>' " • (3) 

k 

The real and positive parameters are called decoherence parameters and determine the strength of 
the interaction. 

In the following treatments we work with a two qubit system where the Hilbert space is given hy Ti — 
(8) C^. We require the Lindblad generators Pk to project onto one-dimensional subspaces and satisfy 
"^k^i Pk — 1- Furthermore we constrain to one dissipation parameter A which parameterizes the strength 
of the interaction and therefore of the decoherence [l^. Then the dissipator can be written in a very 
compact way 

4 

V{p) = \{p-Y,PkpPk) , (4) 
which is easier to deal with than the original one. 



B. Spin geometry 

The state space for a system consisting of two qubits has in general four complex or eight real dimensions. 
A general density matrix of such a system can be expressed in the following way Q 

p — i(l ®l + rh(j®l + nl®a + Cij aj) , (5) 

where 0-^(8)1, l^aj, cTii^aj} forms a basis in terms of product Pauli operators of S(7i) = S(C^®C^), 

the algebra of bounded operators on the Hilbert space H.. 

The local parameters m, n e R'^ of the density matrix ([5]), given by to Tr {(f(^ 1 p) and n = Tr {1(^(7 p), 
determine the reduced density matrices, e.g., pi = Tv 2P = ^(1 + to • a), whereas the real 3x3 matrix 
c = (cy), where Cy = Tr (di ® aj p), determines the nonlocal correlations. The expectation value for a 
joint spin measurement in directions a and /?, given by E{d, 13) = Tr [p a ■ a ® (3 • a) — {d,c j3), is fully 
determined by the correlation matrix c. 

A state is called separable if it can be written as a sum over product states, that means Cy = • rrij. 
Nonseparable states are called entangled (see, e.g., [13 )■ The amount of entanglement is invariant under 
unitary transformations of the form Ui ®U2- This provides us with an equivalence relation for states 
with equal properties concerning separability and entanglement and we only have to choose a proper 
representative to reduce the number of parameters needed to describe the states. 
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The parameters to, n and c of p transform under the action of the unitary transformation Ui <S) U2 as 
to' — Oifn, n' — and c' = OicO|". The transformation Oi (O2) is related to Ui {U2) via the 
homomorphism connecting the groups SU(2) and SO(3): for every unitary transformation U G SU{2) 
there exists a unique rotation O e 50 (3) such that [/n • all'' — {On) ■ a. 

It turns out that the orthogonal transformations can be chosen such that the matrix c' is diagonal. 
Thus it is sufficient to consider states as representatives where the c-matrix is diagonal (singular value 
decomposition). The singular values, which are different from the eigenvalues, are always real ^. They 
can be arranged to form a 3 dimensional vector c = (ci, 02,03)^ which we call correlation vector. The 
spin geometry picture consists of all possible correlation vectors c. 

Example. 

The correlation vectors for the 4 maximally entangled Bell states, l^*^) ~ '^^^^^'^ 1-'^'^)) "^ith c — 
(±1,±1,-1)^ and = 7f(|00) ± 111)) with c = (±1, =f1, +1)'^, form the corners of a tetrahedron 

which characterizes the convex set of all possible states. The partial transposition condition leads 
to a reflection of the tetrahedron and the set of separable state is given by the intersection of both 
tetrahedrons which results in an octahedron (see Fig. [1] and Ref. jll). 



1^ 




FIG. 1: The tetrahedron of possible states, the inverted tetrahedron and the octahedron of separable states. 

State with the same purity, measured hy S — Tr (p^), correspond to a sphere of a radius proportional to 
i5 centered at the origin of the tetrahedron. There are 4 pure states with 5 = 1 in the picture, the Bell 
states. States with equal entanglement, measured by the concurrence C [l3l |. are found on planes with 
normal vector equal to (—1, —1, —1) for instance. The only maximally entangled states in this picture 
are the Bell states. The separable states (C = 0) lie at the border of the octahedron and inside. For a 
detailed analysis of these facts see Ref. IJ] . 



C. Singular value decomposition 

The singular value decomposition (SVD) is a generalization of the eigenvalue decomposition which can 
be used to diagonalize rectangular matrices (the eigen decomposition is defined only for square matrices). 
In analogy with the eigen decomposition, which decomposes a matrix into two simple matrices, the main 
idea of the SVD is to decompose a rectangular matrix A into three simple matrices: two orthogonal 
matrices P and Q and one diagonal matrix D, given by 

A^PDQ^ , (6) 

where P^P = 1 and Q'^Q — 1. The diagonal entries of D are called singular values. 



Note that pure unitary dynamics do not affect the degree of entanglement and thus do not appear in the singular value 
decomposition nor in the spin geometry picture. For this reason it is legitimate to consider only the pure decoherence 
part of the equation, see Eq. (0. 
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From the relations AA^ = P P^ and A^ A — Q it is obvious that that the eigenvalues of 

the matrices AA^ and A^ A correspond to the squares of the singular values of A and the eigenstates of 
AA^ and A form the columns of the matrices P and Q. The columns of P {Q) are called left (right) 
singular vectors of A. 



II. DECOHERENCE MODES 



The following investigations are concerned with decoherence effects. The aim is to get information on 
the time evolution of the parameters to, n and c = (cy), introduced in Eq. (O, and to construct the spin 
geometry picture. Because the correlation vectors consists of the singular values of the correlation matrix 
where the dynamical evolution has no impact on we can base our investigations on the master equation 

|p(<) = -V{p(t)) = \{p~j2 P^PP^) ' (7) 

k=l 

which describes only effects arising due to pure decoherence and includes no dynamical effects. 

For a two qubit system {|ei;)}fc=i 4 denotes an eigenbasis of the Hamiltonian of the undisturbed system, 

H\ek) = -Efclefc). In the following we chose several kinds of projection operators as Lindblad generators 
for the dissipator, Eq. ([7]). 

A one-dimcnsional (bilocal) projection operator Pk on the eigenspace of the Hamiltonian H can be written 
in terms of projections on the individual subspaces, e.g., 

- 1.= , I - p(l) ^ p(2) 



P^^\eu){eu\^Pl,''®Pl'\ (8) 



where Pi^^ — { ' " " and P^f' = \ ^ , 

[i(l-cr^) for A: = 3, 4 ^ |i(l-cr^) forfc = 2,4 

Decoherence modes where the Lindblad generators are projections of the form ^ are denoted as mode 
A. We consider also (bi-)unitary rotations of the projection operators of the form 

Pk = (Ui ® l)PkiUi ® 1)^ = UiPI'-^uI ® Pf ^ , (9) 
which are called mode B, and more generally 

Pk = {Ui ^ U2)Pk{Ui ® U2)^ = UxP\^^ul ® C/aPf , (10) 

which denote mode C. It is clear that mode A and B are just special cases of mode C. 
A unitary rotation U G SU(2) can be written in the following way 

a a / cos § — ia^ sin ^ — (m^; + a„) sin ^\ 

C/ = e-'^'^-" =cos^l-«sin^a-a= ' ' , (11) 

2 2 \ _ Q,^) sin ^ cos ^ + sin ^ / 

where a denotes the rotation angle and the unit vector a indicates the axis of rotation. For convenience 
we chose a = (0, 1, 0)^. Thus the rotation matrix has the simple form 



Ul + a,) forA: = l,2 , 0(2) _ j ^(1 + forfc = l,3 



Vsinf 



cos ■2 



and we set Ui = U{a) and U2 = C/(/3) in Eqs. Q and (fTUll . 



A. Mode A 



The first mode describes the simplest possible case. We choose the Lindblad generators according to Eq. 
((HI). We have to calculate the expression Pj^pPj^ which gives the following, based on Eq. 

P^pP^ = \{p^ + rf^ {Pi'^<^Pi'^)®P^k'+nPi'^®{Pl^^SPl^^) + c.,, {Pi'^a,P^'^)^iPl:'^<y,Ph^ ■ (13) 
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After a short calculation we find that 

p(l)^. p(l) _ ^ , ^fc=l,2 g p(l) p(2)„.p(2) _ ^ , ^fc=1.3 g p{2) 

where the signs are chosen accordingly and 5ij denotes the Kronecker-delta, i.e., 5ij = 1 for i ^ j otherwise 
it is 0. The sum over all k gives 

Y,P^pP^ = -(l + m,, {P^ + P^-Pi-Pt)+n,. {P^-P^ + Pi-P^)+c.,. {P,^-P^-Pi + Pt) 
k ^ 



(15) 



= i ( 1 + (g) 1 + 1 (g) + CTz ® cTz ) , 

where we have used the identities 

Pi + Pi - P3 -Pi=(T,®i, 

pA _ pA ^ pA _ = 1 ^ , (16) 
P^ - P^ - P^ + P^ = G,®G,. 

With Eqs. (HI and (fT^ we are ready to evaluate the master equation A comparison of coefficients 
leads to the following differential equations 

\ / \ / \ / '^x \ / ^xx ^xy ^xz \ / (^xx ^xy ^xz \ 

I = — A I I , I I = -A I Uj^ I , \CyxCyy Cyz l = ~ l ^^V^ Cyy Cy z I , ( f 7) 

^lilz J V / J \0 J \Czx Czy Czz J \Czx Czy / 

with the solutions 

mz{t) = m,(0), nz{t) = n,(0), Czz{t) = €^^(0) 

mj(t) = e"^*mj(0), n,(t) ^ er^^n,{0) fori 7^ z (18) 
c,j{t) = e~^*c,y (0) fori = j 7^ z . 

This means that all elements are damped by e~^* except the z-components of the local parameters and 
the zz-coniponent of the correlation matrix which are unaltered. 

B. Mode B 

The next stage is to consider rotations in one subspace, e.g., projection operators of the form Now 
the expression P^ pP^ looks like 



Pk''pPi = \{Pk^+r^ {UiPl^M<^U,Pl^M)®P^k^+n (^iPf^t^l)®(Pf a/'f ^ 



(19) 



With the expression 

UiPI^^uIg.UiPI^^uI = (±)fcZ3;4 (<5„ sina + 8,^ cos a)UiP^'^ Ul , (20) 

we get for the sum 

P^ pP^ = i ( 1 + {rrix sin a + iriz cos a) (sin acr^ 1 + cos a CTz 1) + ?iz 1 (g) cr^ 



{cxz sin a + Czz cos a) (sin a az + cos aaz Tz) 



(21) 
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The differential equations for the parameters of the density matrix are given by 

?7ir — rrix sin^ a — niz sin a cos a 
my 

— TUx sin a cos a — cos^ a 









rhy 


]- 


-'{ 


\mz^ 










^xy 


I 


Cyx 


'^yy 






Czy 





( rix^ 




]■ 




-A 









/c 



■XX ^xy ^xz 

Cyx Cyy 

Czx Czy C-zz 



Cxz sm a - 

Cy 



Czz sma cosa> 



Cxz sm a cos a — Czz cos a , 



(22) 



Apart from the solutions already given in Eq. (|T8|) we get the following solutions for the remaining 
components 



mx{t) ^ {sm a + e cos 0)7712; (0) + (1 — e ) sina cosa TO2(0) , 
mz{t) = (1 — e"'*'*) sinacosaTO2:(0) + (e"'*'* sin^ a + cos^ a)mz{0) , 



(23) 



Cxzit) = {sin^ a + e cos^ a)c2:z(0) + (1 — e '^*) sin a cos a 0^2(0) , 
Czz(i) = (1 — e^'^*) sin a cos a Ca;z (0) + (e^'^* sin^ a + cos^ a)cz2(0) . 



(24) 



This means the rotation of the projection operators results in a coupling of different components, in this 

case between rux and and Cxz and Czz- 

This mode has been discussed in detail in Ref . [121 . 



C. Mode C 

Independent rotations in both subspaces, where the projection operators are given by Eq. pop . results 
in an expression like 



1 



-C, (C/iPi'^C/iV.C/iPi^V/)®(t/2Pf ;72V,C/2Pf c/t)) , 



PkPPk = 4 ( Pk+^ {U^Pl'>UlaU,Pl'>Ul) ® {U^P^Ul) + n {U.Pl^'Ul) ® {U2Pl^^UlaU2PrUl) 



3(2)7-rt^ 



(25) 



(26) 



where UiP^^'^Ula.UiPl^^Ul is given in Eq. ^ and 

C/2Pf' C/la.C/aPf = (i)^ {S^xSm|3 + 5,z cos/3)[/2Pf . 
Calculating the sum over all 4 terms gives 

PkpPk = ^ f 1 + m (sin a (Tj; (g) 1 + cos a CTz (g) 1) + n (sin (3 1 (g) ax + cos /3 1 (g) az) 
k ^ 

+ c (sin a sin (3 ax ® <7x + sin a cos /3 CTx ^ Cz + cos a sin (3 az ® ax + cos a cos I3az® az 

(27) 



where 



m = nix sin a + Wz cos a , n = Ux sin P + Uz cos /3 , 

c = Ca;a; sin a sin /3 + Cxz sin a cos /3 + Cza; cos a sin /3 + Czz cos a cos /3 



(28) 
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The differential equations for tlie parameters of the density matrix are given by 



Thy 



TTLx — rux sm a — niz sm a cos a 
-A I niy 

, TOr — TJix sin a cos a — cos^ a 









1 













nx — nx sin /3 — Uz sin/3 cos/? ' 

Uy 

nz — nx sin j3 cos j3 — Uz cos^ /3 < 



Ca:a: ^2;^ Ca;^ 
Cya; ^yy 



^ Cxx — c sin a sin /3 Ca:j, Ca;2 — c sin a cos (3\ 



^yx 



^yy 



ccosasin/3 Czy Czz — c cos a cos P 



(29) 



The solutions for rn{t) and have the same structure, given by Eq. ([^5]) . The solutions for the modified 
components of the c-matrix are the following 

Cxxit) = e-^*Ca;a;(0) + (1 " e~^*) ( siu^ Q; siu^ ^ c^^a; (0) + sin^ a sin /3 COS ^5 Ca;z (0) 



■ sin a cos a sin (3 Czx (0) + sin a cos a sin /3 cos /3 c^^ (0)1, 



Cxz{t) = e-^'cxz{0) + (1 -e-^*)( sin^asin/?cos/3c:,a;(0) +sin^acos^/3c:,z(0) 

+ sin a cos a sin /? cos /3 Czx (0) + sin a cos a cos^ /3 (0) 
Czx{t) = e~^'cza;(0) + (1 — e~^') ( sin a cos a sin^ /3 Caa;(0) + sin a cos a sin /3 cos /3 Cxz{0) 



(30) 



+ cos a sin /3 Cza;(0) + cos asin/3cos/3 C2z(0) I , 



Czzit) — e Czz{0) + {I ^ e ) sinacosasin/3cos/3 Ca;a;(0) + sinacosacos f3 Cxz{0) 



+ cos^ a sin f3 cos /3 Cza; (0) + cos^ a cos^ /3 c^z (0) 



Now the coupling in the c-matrix is extended to include Ca;a;, Ca;z, Cza; and Czz components. Therefore the 
solutions for this mode are more demanding. 



D. Comments 



We see from the structure of the solutions for mode C, Eq. (jSO)) . that the rotation axis of the unitary 
rotation has an influence on the resulting local parameters and the correlation matrix. The different 
modes correspond to the coupling of the system with different environments. The coupling strength 
parameterized by A is always the same but the effects on the state of the system are different. 
To calculate the correlation vector for all these modes is in general not very easy because one gets quite 
big terms. Therefore we postpone this investigation to the next section where we calculate the vectors 
for the special case of the Bell singlet state. 

An interesting point is what happens with the correlation matrix for t ^ oo. As can be easily checked 
the asymptotic correlation matrix is given by 



sin a sin /? sinacos/3' 
w \ 

,cosasin/3 cos a cos /3y 



(31) 



where w = sina(ca:z(0) cos P + Cxx{0) sin/3) + cos a(czz(0) cos P + Czx{0) sin/3). The correlation vector for 
the asymptotic matrix has the simple structure of Coo = (0, 0, w)-^. This shows that all states independent 
of their initial correlations end up on a line connecting the origin of the spin geometry picture (the 
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maximally mixed state) with a point representing an equal mixture of two Bell states given by the 



corners of the octahedron (see Fig. 3(a) ). Note that the ordering of the singular values does not matter 



due to the high symmetry of the picture. 

III. EXAMPLE: BELL STATE 

To get a better feeling for the decoherence modes we consider the example of the maximally entangled 
Bell singlet state l^*") -^(|e2) — jes)) = -^(0,1,-1,0)^ where the initial density matrix is given 
by p(0) = |4'~)(\E'^|. The local parameters vanish to = n = and the correlation matrix is diagonal 
Cbs(O) = diag(— 1, —1, —1). In the following we consider the behavior of the correlation matrix because 
the local parameters do not change for this state. 

A. Mode A 

For decoherence mode A the correlation matrix remains diagonal but is affected by the decoherence 
like c^g = diag(— e""**', — e""^*, — 1). Consequently the correlation vector for mode A is given by Cgg = 



The correlation vector c^g for fixed A is plotted in Fig 2(a) with respect to varying t. We start in 
the corner indicated by l^*") and for t — > cxd approach the point bisecting the line which connects the 
projectors of the states l^*") and This asymptotic state is mixed but not maximally mixed and 

lies at the border of separability which is given by the blue octahedron. 

B. Mode B 

Mode B, where the projector in one subspace is rotated unitarily, leads to a correlation matrix given by 



-BS 



-e — (1 — e '*'*) sinacosa 

-e--^* I . (32) 

-e-^* - (1 - e"^*)cos2 



The emergence of the a:z-component is due to the particular coupling of the differential equations, Eq. 
The singular values ^ of this matrix are computed to be 



-At j / = (33) 

3 + 2e2^*cos2a - cos(2a) =F v2(e^* - l)cosa^5 - 3cos(2a) + 2e^*(2 



'-2,3 — 

which gives for the correlation vector c^^ — (— cf , — cf , — cf )"^. 

The expansions of the last two singular values up to second order in a, 

-At 1 _ At p-At 1 _ At 

show that to first order in a the singular value C2 is well approximated by e~^* and C3 is constant. The 
second order contributions affect C3 much more than C2 . In fact the deviations of C2 from the exponential 
function are very small (a few percent) and the largest deviation arises for a = 

The decoherence paths for mode B form a plane which has a little bulge due to the slight deviation of 
cf from the exponential function. This deviation from the plane formed by equal cf and cf coefficients 



^ Note, that the SVD does not specify a certain sign, which can be determined, e.g., from the eigenvalue decomposition, 
nor a certain attribution to the coordinate axes, but note the high symmetry of the tetrahedron. 
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can be understood by the fact that the SVD is an asymmetric transformation and although apphed to 
square matrices the matrices P and Q, Eq.®, are different. A comparison with the eigenvalues of Cgg 
given by Af = Af = — e"'*'* and Af = — cos^ a — e 
on the angle a. 

We can distinguish two special cases. The correlation vector for a 
which corresponds exactly to mode A (see Fig 



SVC? a, reveals that only the z-component depends 



c = 



-At 



-At 



-\t\T 



is shown in Fig. 2(b) 



2(a)). For a 



is givenby c = (— e e 1)^ 
: -J the correlation vector, given by 



The state approaches the totally mixed state sitting 
at the origin of the coordinate system. Thereby it reaches the border of separability at At = In 3 (cf. with 
Ref. [l2|)- We recover the phenomenon of "entanglement sudden death" introduced by Yu and Eberly 
[lH. Note that the Werner state [l^, which interpolates between a maximally entangled state and the 
maximally mixed state, shows the same behaviour in the spin geometry picture. 

In Fig. 2(c) the correlation vectors are shown with respect to fixed parameter A and evolving time t for 
The time the border of separability is reached varies with respect to a and the 
= ^ where the border is reached after the shortest time and a = where it is 



different values of a. 
extremal cases are a 
reached asymptotically at infinity. 






(a) 



(b) 



FIG. 2: Tlie red lines represent possible decoherence patlis of mode A and B for the Bell singlet state. Figure 
(a) shows the path traced out by mode A, (b) depicts mode B for Q- = ^ (the hne of the Werner state) and (c) 
visualizes mode B for different a- values in steps of jq. 



C. Mode C 



The most general case, mode C, which is characterized by two bilocal unitary rotations of the Lindblad 
generators, exhibits a time evolved correlation matrix for the Bell singlet state given by 



— e — (1 — e ^*) cos(a — /?) sina sin/3 — (l — e ^*) cos(q! — /?) sinacos/? 

-e-^* 

— (l — e^-^*) cos(a — /3) cosa sin/3 — e""^* — (1 — e^'*'*) cos(a — /3) cos a cos/3 y 



(35) 

We calculate the singular values. 



-1 



e 



-At 



c^3 = ie~^*y3 + 2e2-^*cos2 A - cos(2A) T \/2(e-^* - 1) cos Ay^5 - 3 cos(2A) -I- 2e^*(2 + e^*) cos^ A , 

(36) 

where A = a — /3 and realize that they have the same structure as Eq. (|33p with a symmetry in a and 
(3. They depend only on the relative angle difference |A| and consequently the problem reduces to mode 
B, already discussed in section UlIBI 



Just for completeness, the eigenvalues of the correlation matrix for mode C are given by Af = A2' = 



-e 



-At 



and Af = - cos^ A - e~^* sin^ A. 
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IV. EQUIVALENCE OF MODE B AND MODE C IN THE GEOMETRIC PICTURE 

The results of the last section suggest to investigate the point which initial states result in an equiv- 
alence of mode B and C in the spin geometry picture. This means that for some initial conditions 
the singular values (and also the eigenvalues) of the correlations matrix for mode C depend only 
on the angle difference A = a — /? and have the same structure as those for mode B. We state 
a proposition which is a necessary condition for the initial correlation matrix that mode C depends 
only on the angle difference and is equivalent to mode B. We do not know if the condition is also sufficient. 

Proposition. 

For a density matrix with initial local parameters m(0) and n(0) equal to zero and an initial correlation 
matrix given by 

/ fcl k2\ 

c(0) = ^3 , (37) 
\-k2 kj 

the geometric picture of decoherence mode B and mode C coincide and the singular val- 
ues depend only on the angle difference A. The correlation vector of matrix (|37p is given by 

c(0) ~ {\/kl + /c| , \Jk\ + k2 , fca) where the values fci, k2 an d k^ h ave to satisfy 2-y/fcf + fcf + fcs < 1 
in order to belong to the tetrahedron of possible states (see Fig. |3(b)[). 



Proof. 

We consider states where the initial local parameters are set to zero because of their irrelevance for the 
spin geometry picture. 

Two matrices A and B have the same eigenvalues Xi if they are similar A = SBS^^ for an invertible 
matrix S. That means two conditions have to be satisfied: 

• TtA = TtB = J2^^^ , 

• dot A = det B = n,: ■ 

In our case we have to deal with the matrices A"'" A and B^ B and find out when their eigenvalues are 
equal. In general the equality of the singular values of two matrices does not imply the equality of the 
eigenvalues of these matrices, e.g., in general 

A^A = (SBS-^fiSBS-^) ^ S{B'^B)S-^ . (38) 

The relation A^A = S{B^ B)S^^ is only valid for orthogonal transformation matrices S where S^^ — S"^ . 
Let us assume that this condition is valid, then it is sufficient to consider only the conditions stated above^. 
In the following we set A and B c*^. 

There is only one point left which is important. The functions of A depend only on the angle a which 
we have to substitute with the difference a — f3 in order to compare it with the functions of B. Then we 
calculate the difference of both solutions and determine the conditions for it to be zero and independent 
of the rotation angles. 
We get for the trace condition 

TrB - Tr A = (1 - e"^*) sin/3 cos(a - l3)(^{cxz + c^x) cosa + {c^x " c^^) sina^ . (39) 

This product is zero and independent of the single rotation angles only for 

Cxx = = fcl and c^z = -Czx = ^2 • (40) 



^ We have checked that Tr (A^A) and det(A^A), in contrast to TrA and det A, does not reveal new constraints on the 
coefficients. 
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All other entries of the correlation matrix are not fixed by this condition and can be chosen arbitrarily. 
The condition of the determinant is given by 



det-B - detA =e"^*^ (l - e"^*) smp(ki cos(a - /3) + fcg sin(a - /3)j 

(j^^xyCyz ~^ Cyx^zy) COS O. -f- {CxyCyx ^yz^zy') OL^ , 



(41) 



where we already used the condition stated in . This is equal to zero and independent of the rotation 
angles for 







(42) 



Thus we have only one entry, namely , in the initial correlation matrix left which is not determined by 
both conditions. □ 

Note that the invariants introduced in Refs.[l3, [iBl are connected with the problem considered here. 
The states given by Eq. (|37p form planes (with proper restrictions of the coordinates) which contain 
the lines connecting the Bell states and the line of the Werner state (see Fig. 2(b) I in the spin geometry 
picture. 





(a) 



(b) 



FIG. 3: Figure (a) shows the asymptotic states of mode C in the spin geometry picture, given by Eq. (|3ip . In 
|(b)| all states which result in an equivalence of mode B and C are depicted. Note the symmetry of the picture of 
which we have picked out only one possibility. 



SUMMARY AND CONCLUSION 



We considered the master equation with a special type of dissipator which describes decoherence in a 
two qubit system. The generators of decoherence are chosen to be on the one hand projectors onto the 
eigenstates of the undisturbed Hamiltonian of the system (mode A) and on the other hand they are 
(bi-)local unitary rotations of these projectors (mode B and C). 

The general solutions of the master equation are presented with respect to the decomposition of a two 
qubit system in terms of joint Pauli matrices. We are interested in the time behavior of the correlation 
matrix which is studied for the case of the Bell singlet state. The information about entanglement and 
purity is encoded in the correlation vector which can be illustrated in the spin geometry picture. We 
show graphically the paths of the different decoherence modes and discuss their behavior for the Bell 
singlet state (see Fig. (2). 

For the special case of the Bell singlet state we find that mode B and C are equal and depend only on the 
difference of the rotation angles of the projection operators. This arises the question for which general 
initial states this is the case. We conclude that this happens for all states which are contai ned in the 
plane formed by the line connecting the Bell states and the line of the Werner states (see Fig. |3(b)"] ). 
The asymptotic states of decoherence mode C (and therefore also for mode A and B) are found to give 
a line connecting the maximally mixed state at the origin of the picture with the equal mixture of two 
Bell states (see Fig. 3(a) I. 
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The decoherence modes investigated in this paper exhibit very strong symmetry properties which is due 
to the fact that we choose projection operators as generators of the decoherence. For more general 
Lindblad operators the calculations start to get more involved but the conjecture is that for them the 
whole geometric state space of the tetrahedron can be occupied. 

An open question related with the spin geometric picture is the influence of the local parameters on the 
geometric state space? We have some preliminary results for maximally entangled mixed states [l^, [l^l 
where the local parameters are not zero any more and the calculation of the singular values is more 
involved. 
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